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The glueball-to- vacuum matrix elements of local gluonic operators in scalar, tensor, and pseudoscalar channels 
are investigated numerically on several anisotropic lattices with the spatial lattice spacing in the range 0.1 fm 
- 0.2 fm. These matrix elements are needed to predict the glueball branching ratios in J/tjj radiative decays 
which will help to identify the glueball states in experiments. Two types of improved local gluonic operators 
are constructed for a self-consistent check, and the finite volume effects are also studied. The lattice spacing 
dependence of our results is very small and the continuum limits are reliably extrapolated. 



1. Introduction 

Glueballs, predicted by QCD, are so exotic 
from the point of view of the naive quark model 
that their existence will be a confirmation of 
QCD. Extensive numerical studies have been car- 
ried out to simulate the glueball spectrum and 
resulted that the low-lying glueballs are in the 
mass range 1-3 GeV, which suggests that the 
J I ip radiative decays are the ideal hunting ground 
for glueballs. There are several possible glue- 
ball candidates in the final states of Jji\) radia- 
tive decays, however, more criteria are needed for 
their unambiguous identifications, one of which 
might be the partial widths of J/ip radiative de- 
cays into glueballs. To estimate these partial 
widths, the vacuum-to-glueball transition matrix 
elements (TME) of local gluonic operators should 



*Presented by S.-J. Dong. 



be derived first. 

The techniques of lattice simulations in the 
glueball sector have been substantially improved 
in the past decade. Inspired by the success of 
anisotropic lattice techniques in the simulations 
of the glueball spectrum pQ and as a continua- 
tion of former studies |2I3| , this work is devoted 
to the numerical study of TME on anisotropic 
lattices with tadpole-improved gauge action. 

2. LOCAL GLUONIC OPERATORS 

The TME computed in this work are 
(0|S(aO|0++), (0|T MV (x)|2++), and (0|P(x)|0-+), 
where \J PC ) refers to the glueball state with 
the quantum number J PC , and the local oper- 
ators S(x), T^, and P{x) are trace anomaly 
g 2 TrG p „G llv (x), the energy-momentum tensor 
g 2 Tr(G l _ ia {x)G al ,(x) - \g ^G 2 (x)) , and the topo- 
logical charge density g 2 € pypcr TrG \ LU {x)G pa (x), 
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Table 1 



The simulation parameters. 









a s (fm) 


L A xT 


#Meas. 


2.4 


5 


0.409 


0.222(2) 


8 a x 40 


20000 










12 3 x 64 


10000 










16 3 x 80 


10000 


2.6 


■5 


0.438 


0.176(1) 


12 3 x 64 


8600 


2.7 


5 


0.451 


0.156(1) 


12 3 x 64 


10000 


3.0 


3 


0.500 


0.120(1) 


16 3 x 48 


10000 


3.2 


3 


0.521 


0.101(1) 


24 3 x 72 


7900 



respectively (G M „(x) is the gauge field strength 
and g the gauge coupling). Two types of lattice 
local gluonic operators (Type-I and Type-II) are 
constructed in this work. 

Type-I operators are linear combinations of a 
set of small Wilson loops according to the irre- 
ducible representations of the lattice symmetry 
group, namely, Af + , Aj + , E ++ , and T 2 ++ gj. 
Letting O r (x) be the local operators of specific 
quantum number R, the construction of Type-I 
operators can be expressed as 

R {x) = C*ReTr(W± l \x)+aW¥\x)+. . .), (1) 

i 

where Cf- are the combinational coefficients. Dif- 
ferent Wilson loops W\ , W2 , ■ • ■ are included with 
proper factors a to improve the operator. 

To construct Type-II operators, we define the 
lattice gauge field strength F flu {x) as 

F^{x) = Im(f(u s )P^{x) + g{u s )R^{x)) Cl (2) 

where P fll/ (x), and R^(x) are respectively the 
plaquette and rectangle at x. (|) c means the 
clover average, f(u s ) and g(u s ) are factors includ- 
ing the tadpole parameter u s — ((l/STrPij)) 1 / 4 
and are chosen in a way that yields F^ v {x) — 
a l(G^{x) + 0(a 4 )) (a s is the spatial lattice spac- 
ing). Proper combinations of TrF fll/ (x)F pa (x) 
give the Type-II operators with reduced lattice 
artifacts. 

There are seven matrix elements calculated in 
this work, which are denoted by (S,B), (S,E), 
(E,B), (E,E), (T,B), (T,E), and (PS). Here ( ,B) 
(or ( ,E)) means the operator is made up of color- 
magnetic (or color-electric) field, and (PS) refers 
to the pseudoscalar channel. 
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Figure 1. Finite- volume effects of matrix ele- 
ments. 



3. SIMULATION DETAILS AND RE- 
SULTS 

Numerical simulations were carried out on 
anisotropic lattices with tadpole-improved gauge 
action £Q. Five independent simulations have 
been done with input parameter listed in Table. 
1. 

With the tadpole improvement, the local glu- 
onic operators on the lattice are all improved to 
0(a*) (a s is the spatial lattice spacing) at the tree 
level. With applying the variational method to 
the combinations of smeared Wilson loops with 
different prototypes P, glueball states are ob- 
tained through correlators of smeared operators 
which have large overlaps with glueball states. 
Six independent runs were carried out, on lat- 
tices with spatial lattice spacings in the range 0.1 
fm — 0.22 /m, to measure the smeared-smeared 
correlators Css(t) and smeared-local correlators 
CsL{t). The matrix elements are extracted by fit- 
ting Css(t) and Csh{t) simultaneously using the 
correlated x 2 method. The fit models are taken 
as 

C SS (t) = X 2 e- Mt (3) 
C SL (t) = XYe~ Mt , 

where X is the amplitude of the glueball opera- 
tors, Y is the glueball-to- vacuum matrix element, 
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Figure 2. Continuum extrapolations of matrix 
elements of Type-I operators. 
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Figure 3. Continuum extrapolations of matrix 
elements of Type-II operators. 



and M is the glueball mass. 

The finite volume effects (FVE) of matrix ele- 
ments are studied on lattices 8 3 x 40, 12 3 x 64, 
and 16 3 x 80 at (3 = 2.4, f = 5, and the results are 
shown in Fig. ^ where each point with a errorbar 
is the fractional change, Sg(L) = 1 — fc(L)/ /g, 
in the matrix elements (fa is the average value of 
the matrix elements of glueball G over the three 
lattice volumes, and fc(L) is the matrix element 
of glueball G measured on lattice L 3 x T). The 
labels L = 8, 12, and 16 denote the different lat- 
tice volumes, and the labels along the vertical 
axis represent the matrix elements of the differ- 
ent local operators. To guide eyes, Sq = and 
Sq = ±0.02 are also drawn in Fig. ^with solid 
line and dash lines, respectively. All changes are 
statistically consistent with zero, indicating that 
systematic errors from FVE are negligible. 

The simulated results of matrix elements at dif- 
ferent lattice spacings are shown in Fig. [3 (Type-I 
operators) and Fig. (Type-II operators), where 
a s (in units of r ~ (410(20) M eV)' 1 , the 
hadronic scale parameter) dependences are ob- 
served. The matrix elements can be reliably ex- 
trapolated to the continuum limit by using the 
form f{a s /r ) = Trjj + c 2 (a s /r ) 2 + c 4 (a s /r ) 4 , 
where TYg, C2, and C4 are best-fit parameters. 
We keep the a 2 s term in the fitting model because 
there are residual ot 8 a 2 s artifacts in the gauge ac- 



tion and some local operators. From the figures 
one can find that the Type-II operators exhibit 
better behaviors (for example, the matrix ele- 
ments of T2 and E representations coincide, as 
it should be when the rotational symmetry is re- 
stored in the continuum limit). 

The physically available predictions will not be 
derived until the local gluonic operators are prop- 
erly renormalized. The nonperturbative renor- 
malization of these operators is in progress. 
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